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La presente invention conceme une methode et un dispositif de detection multi- 
utilisateur. Plus particulierement, la presente invention conceme une methode et un 
dispositif de detection multi-utilisateur pour un systeme de telecommunication DS- 

CDMA ou MC-CDMA. 

5 Dans un systeme de telecommunication mobile DS-CDMA (Direct Sequence 

Code Division Multiple Access), la separation des communications en provenance ou k 
destination des diflferents utiUsateurs est realisee en multipliant chaque symbole 
complexe d'un utilisateur par une sequence d'etalement qui est propre h ce dernier, 
encore appelee pour cette raison signature de I'utilisateur. La frequence d'etalement 

10 (chip rate) etant sup6rieure a la frequence des symboles, le signal transmis par chaque 
utilisateur est distribue (ou 6tale) dans I'espace des frequences. Le rapport entre la 
bande occupee par le signal etale et la bande occup6e par le signal d'information est 
appele fecteur d'etalement. A la reception, la separation d'un utUisateur donn^ est 
obtenue grace a un filtrage adapte a la signature correspondante. Lorsque le canal de 

15 transmission presente une pluraUte de trajets de propagation, la sortie du filtrage 
adapte comporte autant de pics de correlation. Chaque trajet du canal peut etre 
modelis6 par un coefficient multiplicatif complexe et un retard. Les signaux s' etant 
propages selon les dififerents trajets peuvent etre alignes et combines au moyen de 
coefficients complexes conjugues des coefficients de trajet, realisant ainsi un filtrage 

20 adapte au canal de transmission. Pour simplifier la terminologie nous engloberons 
dans I'expression generale « filtrage adapte a I'utilisateur k » h ]a. fois I'operation de 
filtrage adapte a la signature de I'utilisateur k et celle de filtrage adapts au canal de 
transmission. 

Pour combattre I'interference entre signaux h destination (lidson descendante) 
25 ou en provenance (liaison montante) des diflferents utilisateurs, il a 6t6 propose des 
methodes de detection multi-utilisateur et notamment des methodes de detection 
itdratives telles que celles connues sous le nom de PIC (Parallel Interference 
Cancellation) et SIC (^Serial Interference Cancellation). Elles sont fondees sur 
Titration d'un cycle d'elimination des interferences coii5)ortant 1' estimation des 
30 symboles emis, revaluation des kiterferences et leur soustraction des signaux re9us. 
Bien que performantes, ces methodes ne sont pas optimales dans la mesure ou elles ne 
foumissent pas une estimation au sens du maximum de vraisemblance des symboles 
transmis par les differents utilisateurs. 
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Une methode de detection multi-utilisateur a maximum de vraisemblance 
inspiree de I'algorithme de Viterbi a ete proposee par S. Verdu dans un article intitule 
« Minimum probability of error for asynchronous Gaussian multiple access channels » 
publie dans IEEE Transactions on Information Theory, pages 85-96, Janvier 1986 
mais sa complexite est prohibitive car eUe varie de maniere exponentieUe avec le 
nombre d'utilisateurs . 

Plus recemment il a ete propose par L. Brunei et al., dans un article intitule 
« Euclidian space lattice decoding for joint detection in CDMA system » public dans 
Proceedings of ITW. page 129, Juin 1999. une methode de detection multi-utilisateur 
a maximum de vraisembbnce utilisant une representation par reseau de points. Selon 
cette methode on determine un vecteur caracteristique du signal re9u representant une 
statistique sufBsante pour la detection au maximum de vraisemblance des symboles 
transmis par les differents utilisateurs. On montre sous certaines conditions que le 
vecteur caracteristique pent etre repr&ente comme le point d'un reseau perturbe par 
un-bruit. La detection consiste alors-a rechercher le-point-du-re 



point correspondant au vecteur re9u. Cependant, la dimension du reseau k utiliser 6tant 
en general de 2.K ou K est le nombre d'utilisateurs le nombre de points a tester est 
eleve. Pour simplifier la detection, U a ete propose de limiter la recherche du phjs 
proche voisin aux points du reseau appartenant a une sphere centree autour du point 
refu. Nous exposerons ci-apres cette methode de detection simplifi^e dite « methode 
de detection par spheres » : 

Nous nous pla^ons dans le contexte d'un systeme de telecommunication mobUe 
k acces multiple avec 6talement de spectre par sequences directes (DS-CDMA) 
comprenant K utilisateurs communiquant de maniere synchrone avec une station de 
base. 

Soit diiO le symbole complexe ^mis par I'utilisateur k a I'instant /. Ce symbole 
appartient k la consteDation de modulation Ak utilis^e par I'utilisateur k, que I'on 
appellera encore alphabet de symboles de Tutilisateur k. Cheque utilisateur k 
transmet un bloc de N symboles avec une ampUtude du signal m. Les symboles sont 
etales par une signature complexe s,(t) = sUt) + j\sUt) de dur^ egale a la p^riode 
symbole T : 

s,(O = 0siti[0,T[ 




Les K symboles complexes dkii)=dk{{)+j.dk{i)\xBi)STcas a I'instant / sont places dans 
vin vecteur ligne de valeirrs reelles d2(i) defini comme : 

dii})={d^{i)4i{j\ . ^d^(i),dk(i)) (1) 

Le signal module correspondiant est alors, en fonction du temps / : 

S.JyyMd,ii}sk(t-iT) (2) 



/«0 k=} 



Nous supposons que le canal est un canal ideal a bruit blanc additif gaussien. Soit 
n=St+7jf le signal re9U au temps / et 7, un bruit gaussien complexe de moyenne nulle 
et dont les composantes ont une variance . 

Soit le vecteur Hgne y2(f)4}^i^(f^y{(f)^^^y^^^ tel que j*(0=yKO+y>/(0 soit la 
10 sortie complexe a Tinstant / du filtre adapte a Tutilisateur k : 

yt(i)= ^k\t-iT)ndt 

—00 

=Ja.J/(Ot/(/>*(/)d/-fm(0 (3) 

/-I 0 

=^aidi{i)Rik-^m{i) 

T ^ 

avec Rtk^Uiit)sk\t)6t=Rik+jJiik pour A, ^ = 1 K ct m{i)=y^k\t-uT)^ 

La matrice R d' autocorrelation des sequences d*etalement depend en general de 
rindice i car les signatures des diflFerents utiUsateurs peuvent changer au cours du 
15 temps. Cependant, pour des raisons de simplification de notation, nous omettrons cet 
indice par la suite. 

Si Ton decompose les elements complexes de (3) en leurs parties reelles et 
imaginaires, on obtient : 

20 [y«/)+ j.ylOihfl^'ibt^^^^ 



Soient A2=Diag(ai,ai^.^Jc,aif) et R2 la matrice de taille 2K^2K telle que : 



R2= 



-RkxMf—RkKRh^ 



(5) 
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L-aquation (4) peut alors se mettre sous forme matricielle 
y2(i>d2(/)M2+n2(/) 



(6) 



10 



' "VM.": VWVJ^ ,nf (iJUiHO) a pour matrice de covariance A',R 
Nous dtooMrerons ci-apres que y,(i) . ,el que do™ par I'equation (6). peut ftre 

veCc^Te 7^ * Poi- A de d™ . .„u. ensemble de 

x=Z»vi+&V2+. . .+&rv, ou AeZ,V/=l , . 
et oii fvi.vz^.^vx:} est une base sur R' . 



20 



25 



Un-exemple de r^seau de points de dimension 2 a 6te represente en Fig 1 
Les points du reseau fonnent un sous-groupe abelien additif de R« c'est 
d'adleurs le plus petit sous-groupe de R' contenant les vecteurs {v,V2..,v.} et un 2. 

^O^pe^™^^ 

a volume fond^nw rfes, autre que le module du produil vectond des ;rveeteu,= 
de tese e, est done ^gal a |de.(G)| od de.(.) d.sig„e le dete^Zt S^TS 

c^Tu^" '"""T ""^ " par 

eontre qu une unique valeur pour le volume fondamental 

r&eau1??r '^JT'^\'^ * d'un point x apparteuant au 

^« «t l««emble des poto de R' phts proches de x que de tout autre point du 

r^eau. U volun« de cette rigion ea <gal au volume fondaiuental 

Le rayon d'«upileme« P du r&eau ea le rayon de la plus grande spM„ 
fflscnte dans la region de Voronol e, le rayon de recouvremem celui de la plus petite 

est eelu. des spheres les plus petites qui. ceMrees sur les points du rfeeau pennettent 
dereeouvrtrtoutrespaee R^Ladensit.dur^ues.,e.rappon.,«,fe.;2::t 
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sphere de rayon p et le volume fondamental. Enfin, le coefficient d*erreur (kissing 
number) r(A) du reseau est le nombre de spheres tangentes a une meme sphere dans 
Tempilement ou, en d'autres termes, le nombre de voisins d'un point du reseau, situes a 
la distance minimale = 2/> . 
5 Considerons a nouveau requation (6). Les composantes du vecteur dzO) 

appartiennent a un alphabet fini A de cardinal : 

Car^(A)=f[CarJ(A*) (8) 

Nous appellerons A la constellation du systeme (ou simplement constellation) par 
opposition aux Ajk dites constellations de modulation. 
10 Supposons par exemple que les composantes et dl{}) soient des symboles de 

modulation PAM d'ordre M : 



^^(/)e{-M+1 . JI/-3 JUf-l} et (9) 

15 J/(/>fM+l -M+3,..,M--3^-l} ' (10) 

Si Ton effectue la transformation : 

rfj«(/)=l(e//'(/)^Af--l) et rfi(/)=^(^i(/)^M~l) soit encore de manidre vectorielle :: 

20 

d'2(/>=^(d2(i>VA/) (11) 



ou VA/=(A/-1^-1,-,A^-1) 

- 25 les composantes et tfi(0 sont des elements de Z et par suite d*2(/) est un 

vecteur deZ^^ 

De maniere generale, s'il existe une transformation affine transformant les 
composantes d^{}) et en des elements de Z, on pent representer le vecteur d'2(/) 

par un vecteur de Z^^. 

30 De maniSre similaire, on eflFectue la transformation correspondante sur y2(/), 

c'est-a-dire : 
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2 V ' (12) 

vecteu'l'X!r' '^^^^t^T'r""!""' ^-VVO..ror. implicite dans la suite, le 

vecteur d.(OM. apparent aloK a un reseau de points Aa de dimension 2 tel que 

5 defini par la relation (7) avec G=M,. Le vecteur y.(/) peut done etre considere 
comme un point du reseau Aa corrompu par un bruit n^C/). 

.le.t t^P'^r <^™P°^tes du vecteur de bruit n.(0 sont des variables 

aleato^es mdependantes gaussiennes centrees. le probJeme de la detection au sens du 

10 rrr"" . ^^"f^^ ^ P- les difKrents utilisateurs se 

lU ramene a la recherche du point du r^«!Mii a t^i 

.ninimale. ^ ^"^^ ^ ^^^^ 

En realite, les composantes du vecteur de bruit n.(/) sont corr^lfes et la matrice 
de covanance de 02(1) est A^^R^ . 

Afinde se i^ner au cas <iecorreMJLFauj^ectuer_^^^^ 

15 une operation deblanchimentdu bruit. uaecoaage 
La Patrice R ^tant hemutienne, la matrice d'autocoirelation R^ est symetrique 
definieposmveet peut done feirel'objetd'une factorisation deCholesky- 

ou W2 est une matrice triangulaire inferieure de taiUe 2J^x2A: ^'^^ 
20 On definitun vecteur d'observationblanchie: y2(/>=y.^)wr' (u^ 

MiiZ Tr^^nr" points constitu. des vecteurs de conpoLes 
t^ l^ct n^ -(/>=x.(/)wr'o. x.(0 est un vecteur de composantes 

appartenant a A.. Le reseau a pour matrice g^n^ratrice 
A2W2, matrice reelle triangulaire inferieure. generatnce 

bruit ffll'^l-^^""' ^^^^^^nt, la matrice de covariance du 

bnnt ffltr. „.(.)Wf est egale aMI. ou I., est la matrice identite de dimension 2K 
detection comprend done une premiere etape de blanchiment du vecteur 

det™ '""^ ^^^'"^'^ ^'"^ ^^^^^ ^" ^" ^^-^ 

0 Pour reduire le nombre de points a tester, comme illustre en Fig. 1 on peat 

luniter la recherche a une sphere centree autour du point y. En pratiqu! le'choix du 
3:1^ . T: ' ' '''' - gld pou" 
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La Fig. 2 represente schematiquement un dispositif de detection multi-utilisateur 
utilisant une methode de detection par spheres. Le signal re9uwest filtre par une 
batterie de filtres adaptes a chacun des utilisateurs, 210i,...,210k. Les composantes 
reeUes et imaginaires du vecteur d'observation y2(/) en sortie des filtres adaptes sont 
5 transmises a une unite de calcul matriciel efFectuant Toperation de blanclument 
spectral selon la relation (14). Les composantes reelles et imaginaires du vectexir 
blanchi y2(0 sont ensuite transmises a une unite de detection par spheres recherchant 
le plus proche voisin du point refu au sein du reseau £l2 de dimension Les 
coordonnees du plus proche voisin donnent directement les composantes reeUes et 

10 imaginaires des symboles estimes dk{i) pour les difKrents utilisateurs. 

L'etape de la recherche du plus proche voisin est gourmande en temps de calcul 
ce qui pent s'averer tres penalisant lorsque le nombre d'utilisateurs est eleve. En outre, 
lorsque le rapport signal a bmit est feible, il est necessaire de choisir un grand rayon 
pour la sphere de recherche et le risque est eleve de devoir considerer un grand 

15 nombre de points du reseau dans Tetape de recherche, 

Le but de la presente invention est de proposer une simplification de la methode 
de detection par spheres notamment en presence de faible rapport signal a bruit. 

L'invention est definie par une methode de detection d'une pluralite K de 
symboles (dk(0) transmis par ou pour une pluralite K d' utilisateurs a partir d'un 

20 signal re9u, chaque symbole d'un utilisateur appartenant a une constellation de 
modulation, la methode de detection mettant en oeuvre un reseau de points (H) 
genere par lesdites constellations de modulation, ladite plurality de symboles des 
difFerents utilisateurs etant representee par un point parmi un sous-ensemble de 
points dudit reseau, dit constellation et le signal re9u etant represente par un point 

25 caracteristique de ce signal, dit point re9U, translate d'un point de ladite constellation 
par im vecteur de bruit (n), la methode comprenant une etape de projection 
orthogonale du point re9U sur vin sous-espace a£5ne, dit de projection, parallele ou 
confondu avec un sous-espace a£Qne delimitant ladite constellation et ime etape de 
recherche du plus proche voisin du point ainsi projete parmi les points de ladite 

30 constellation. 

Avantageusement, ledit sous-espace affine delimitant la constellation est 
determine en fonction de la position du point re9u par rapport a ladite constellation. 
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La recherche du plus proche voisin est avantageusement Mr. J 

par une elape de ffltrage adaple du sigral iwu le filta„. 
—on e. aux .igna... des dmS:!.::,;^,:!:^^ 

La Fig. 1 represente un rtseau de pomts utile 4 b „^i, ^ . 
«nploy&dansler«cep.e„rilh.sWenFig 2- ° ^ '^«'thn 

I- Fig. 5 ^ ^ ^ ^^^^ ^ ^^^^^^ 
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Considerons a nouveau un systeme de teleconununication DS-CDMA avec K 
utilisateurs fonctionnant en mode synchrone. Comme on Ta vu plus haut, la detection 
des symboles transmis par les differents utilisateurs au sens du maximum de 
vraisemblance pent se ramener a la recherche parmi les points d'un reseau (Qa) du 
5 plus proche voisin du point correspondant au signal rcfu. 

Dans le cas ou les sequences d'etalement sont reelles ou plus generalement 
multiples reelles d'un mSme nombre complexe, on pent montrer que la recherche peut 
etre effectuee dans un r6seau de dimension reduite k IC effet, les termes 
imaginaires de la matrice R2 et par suite de la matrice M2 etant nuls, on peut se 
10 ramener a un reseau de points reel A de dimension K et de matrice generatrice M : 



15 ou y^(/),d^(/),n^(/) (resp. y^(0,d'(i),n'(/) ) sont les vecteurs const itues des parties 
reelles (resp- des parties imaginaires) des composantes de y(0,d(0,ii(/) ; 



M=ARou R est la matrice constituee par les coefficients i?/*= k/(/>*(/)d/ et A est le 



vecteur des amplitudes des K utilisateurs. Les vecteurs d'observation y^{i) et y'(0 

appartiennent a R^. Apres transformation eventuelle selon une relation du type de 
20 celle de (12), les vecteurs y^(/) et y^(Opeuvent etre consideres comme des points d'vm 

reseau A de matrice generatrice M corrompus par du bruit. 

On montre aisement que les vecteurs de bruit n^(/) et n' (/) ont tous deux pour 
matrice de covariance NqSSJ) . R etant une matrice symetrique definie positive on 
peut la factoriser selon ime decomposition de Cholesky : R = WW'^ ou W est ime 
25 matrice reelle triangulaire inferieure de taille KxK. ASu de decorreler les composantes 
de bruit, les vecteurs d'observation r6els y\i)^l y'(Osont tout d'abord soumis a une 

operation de blanchiment : 



y'(/)=d'(/)M+n(0 



(15) 
(16) 



T 



0 



30 



(17) 
(18) 



Dans un deuxieme temps, on recherche les plus proches voisins des vecteurs 
y''(i) et y'if) appartenant au r^au de points Q constitu^ des vecteurs 
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fait qu'un symbole d'un utilisateur neut int.rf' T '^'"P*^ 

^^^^ du p^. p.<.,e vois^ au seta <,•„„ r*s=a„^IJr„:^';',rr * 

m^trique ^ ' ""^^'^"^ ^ minimiser la 



metrique 

".ftnque haste sur la matricedecovariance: ^' 1 on uutee 

25 (IS^) 
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L*ensemble des symboles emis est limite a un alphabet de taille finie KKdZr 
appele constellation. Cette constellation est determinee par les constellations de 
modulation utilisees par (ou pour) les a: utilisateurs et le cardinal de Talphabet A^-est 
le produit des cardinaux des difiFerents alphabets de modulation. On supposera que les 
5 points complexes de chacune de ces constellations possedent des vaJeurs reelles et des 
valeurs imaginaires regulierement reparties . 

Comme on I'a vu un decodage exhaustif necessiterait une recherche du plus 
proche voisin dans la totalite de A^r. Le decodeur restreint avantageusement son 
calcul aux points qui sont situes a rinterieur d'une zone de la constellation situee 
10 autoxir du point re9U, preferentiellement a rinterieur d'une sphere de rayon donne 4c 
centree sur le point re9u comme represente en Fig. K Seuls les points du reseau situes 
a une distance quadratique inferieure a C du point re9u sont done consideres pour la 
minimisation de la metrique (19). 

En pratique, le decodeur eflfectue la mirmnisation suivante : 



15 



20 



ni#HK™?JHI (20) 

Pour ce fdre, le decodeur recherche le plus petit vecteur w dans Tensemble 
translate z-3 . On peut exprimer les vecteurs z et w comme : 

z=pG avec p=(/?iv . 

w=^G avec ^=(^,..-y^r) (21) 



II est important de noter que p et ^ sont des vecteurs r6els. Comme w=z-x ou 
25 X appartient au reseau S, on a la relation ^, = p, - P^"^ i=l,...,/c avec w=^4>Vi . 

Le vecteur w est un point du reseau dont les coordonnees 4^,. sont exprimees 
dans le repere translate centre siu" le point re9u. Le vecteur w appartient a ime sphere 
de rayon quadratique C centree en 0 si : 

30 ||wf==Q(^HGG^^"<C (22) 



D2U1S le nouveau systeme de coordonnees defini par ^, la sphere de rayon 
quadratique C centree en y est done transformee en un ellipsoTde centra sur Torigine. 
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La factorisation de Cholesky de la matrice de Gram r=GG'- donne r-AA^ ^ a 
una matrice triangulaire inferieure d'el6ments S>j . ' ^ 

Pre^abla de y, la matnce de Gram T devrait ^tre remplacee par la matrice GR-G^ 
^o"ly~"^ d^finie positive et done susceptible d'une ^torisation dJ 
II faut egalement noter que si le vecteur y a ete blanchi il ti'v « 



15 



Enposanl 



^ (23) 



20 



25 



^»=«S/^ pour /=1 ^ . 

^y=^pour7=l,.,x-^=y+l^.^;j. (24) 
on obtient 

'=> V j (25) 

E„ ^in^res^ .„„ j.^^ 4 la plage de varia«o„s possibles de J, puis e„ 

^ r ^^ ^26) 

«=' I J^i J 

^ monger ,„e les pfi, i^,e« a«x c„„^sa«es e«i^ de . de 
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20 




( K ( K V I 



(27) 



: C- ^ qei ^e-\- Yjl^'^^ +>^+ 
^ t=i^\ \ j=e+\ ) J 



OU \x\ est le plus petit entier superieur au reel xet [xj est le plus grand entier 
inferieur au reel x . 

Le decodeur possede ATCompteurs internes, k savoir un compteur par dimension, 
chaque compteur comptant entre une borne inferieure et superieure comme indique en 
(27), etant entendu qu'a chaque compteur est associe un couple de homes particulier. 
En pratique ces homes peuvent etre mises a jour de maniere recursive. On pose : 

&=S^/^i„..,^.>=p+ ^qji^J (28) 



(29) 
(30) 



T..,=T,^qXS,-h;f 
avec Tk^C 

A I'aide des Equations (28) a (30), on determine r^cursivement la plage de 
variation de chaque composante en commenfant par la composante he : 

n<h<n (31) 

avec Z7=| - | et ZJH J:^^^ | - (32) 



On restreint avantageusement les plages de variation definies par (32) de 
maniere a ne pas tester inutHement des points qui sont situes hors de la constellation. 
On rappelle que chaque utilisateur k utilise une constellation de symholes Aa, chaque 
symhole de la constellation etant un nombre complexe dont les parties reelle et 
imaginaire (eventuellement apres transformation aflSne) sont des elements de Z. 
Pla9ons nous tout d'abord dans le cas general ou le reseau est de dimension 2K^ ou K 
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« ponant le symbote coniplexe (ransmis par ou pour I'milisaleur i r„ 

en Fig. 3. la constellation ou de ™ani*. e,lalel ttlel, T 
modulation de rutiBsateur * e« A;<^ j • " " ""^dbtion de 

i^lle pour la composante conune indique sur la Fig 3 On 

) defimt ensuite un intervalle de recherche [&-.-.,5i-,]par : 

recherche sent ajusties par 0) el 

B -mLM-,).. „^,re . ce ,„e le con,p,eur as^ci. 4 la co„^ f ^ 
P«»ure pas des points qui s„„, situes en dehors de la de rechar h 

-■rrat„To^t-r:r--^^^^^^^^^ 

^ .h^res l«.e » la co„ste«on, Pinter d^The^W^^ 

v»te pour certaines din«„sions i parca que i?<W,- ou bien rlu' ! 

- " 3ph^>;u.rLr - 
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constellation, ce qui peut entrainer de nombreuses etapes de recherche inutiles et done 
considerabiement ralentir Talgorithme. 

Lorsque la constellation a une forme parallelepipedique, par exemple lorsque les 
constellations de modulation des differents utilisateurs sont de type PAM ou QAM, la 
5 presente invention propose de projeter orthogonalement le point re9u sur un sous- 
espace afline (ou variete lineaire) delimitant ou bordant la constellation. La Fig, 4 
represente schematiquement une constellation en forme de parallelogramme (ici /c=2), 
a droite (401) dans la base canonique de Tespace euclidien R*" et a gauche (402) dans 
la base des vecteurs generateurs du reseau. On rappelle qu'un point re9U est represente 
10 par un vecteur z dans la base canonique et par un vecteur p dans la base generatrice, 
tel que z=pG. Dans la representation 401, Tespace est muni de la norme 

euclidienne ||z|f=^2?€5t en 402 de la norme Q(p)=pGG^p^. Ainsi une sphere centree 

sur le point re9u z dans la representation en 401 est transformee en un ellipsoide 
centre sur le point re9u p en 402. Comme illustre en 401, un point re9U n'appartenant 

15 pas a la constellation est d'abord projete orthogonalement sur cette demiere ou plus 
generalement sur une enveloppe E de cette constellation. Le choix du sous-espace 
affine sur lequel Ton effectue la projection depend de la position du point re9u. Ainsi 
par exemple le point Mi appartenant a la zone I delimitee par les droites Di,D2,D4 est 
projete sur la droite Di, le point M2 appartenant au secteur II delimite par les droites 

20 Di et D2 est « projet6 » sur le sommet P12 (c'est-a-dire que Timage de tout point 
appartenant a ce secteur sera le point P12 ), le point M3 appartenant a la zone III 
delimitee par les droites D2,D|,D3 est projet6 sur la droite D2, etc. Comme nous le 
verrons plus loin, ce choix de projection est sous-optimal, II est cependant prefere 
poxir sa facilite d' implementation. Bien entendu, les diflferents types de sous-espaces 

25 affines sur lesquels on peut effectuer une projection dependent de la dimension du 
reseau, Ainsi pour une constellation parallelepipedique generant un reseau de 
dimension 3, les sous-espaces aflBnes a considerer seront les faces, les aretes et les 
sommets du parallelepipede. De maniere generate, pour un reseau de dimension k , les 
sous-espaces aflfines a considerer seront de dimension N=k-L avec \<L<fc, 

30 La projection est plus facile a realiser dans la representation 402. En eflFet, si 

Ton note >9/ les coordonnees d'lm point dans la base des vecteurs generateurs du 

reseau, les droites Di sont definies par les equations : 

>ff/=>ftnin.^=Jl/r-ahiin,r ou bien /Xr=>ftnax,/=M?'+Ofmax,f (34) 



ler depot 



16 



ou et a^, sont des valeurs reelles positives definissant les marges entre la 
— Uon et son enveloppe. On pouna prendre ,es valeurs de ^.et 1. 

efifectuer directement la projection sur la constellation 
' affin. r " 'l-'^'>"q"«> les sous-espaces 

et^e izir r ''"^^ ' "'^'^"^ ^ - 

7 ' """"^'^^ ^ coordonnees fixfes par des 

relations du type (34). On peut supposer sans perte de generality cue T . 
coordo^ees sont les. pren..res et qu'elles sont l^s , la v^^^^e Lf 
10 c'est-a-dire que Pen a p.>j3^, , ^=,,„^ . ""^^^ ' 

Soft z''=p-G la projection sur le sous-espace affine de dimension N du vecteur 

peut done ecnre 

z'=p'G=pV+ltB=zVkB 
IS— o4V=0!U;:i^^.) ea un ve«e»r Hgne^Tdimension i, / ™ vecteu, liglfdr' 
dm««,„ ^. k ™ ™,ew lig« de dimension ^^scntan, les coordonnees iibres 
du sous-espaee veetoriel T associi au sons-espaee afflne 7» B» . T 

(.x«^ con-enan. ies . p.e„«,e, ^ de o'J^ ^ '^^^^ 
cont^n. ,es ^ „ d ^ fbnnan. ™e base du sou^espace yJJr 

20 Autrement dft, on a G = " 



Bj (35) 



done H 'iT'l" « -^-e <ie rang N, on peut 

done definir la matrice pseudo-inverse B^de B On se ran^n. 

vectorieirparsi™p,etranslation.enposant: " """"^^^^ 



25 2''=z''-z'' 

2=z-z° (36) 

On a alors : 



2''=kB=z.B*.B=z.B"(B.B"/'.B 
z''4-zV^.B«)-W 



(37) 



1 er depot 



17 



Le point correspondant au vecteur z est proche de la constellation et sa 
detection sera done beaucoup plus rapide que celle du point initial. Le gain en rapidite 
de detection est surtout appreciable pour les faibles rapports signal a bruit. 

Comme on Ta mentionne plus haut, le choix de projection est sous-optimah Le 
5 choix optimal correspondrait en fait a celui illustre en Fig. 5. Les zones de choix sont 
delimitees par des droites orthogonales aux droites delimitant la constellation. Ainsi, 
les points de la zone I d61imit6e par les droites D'l et D*'i sont projetes sur la droite 
Di, les points du secteur II delimit^ par les droites D'^ et D'2 sont « projetes » sur le 
point P|2, les points de la zone III delimitee par les droites D'2 et D"2 sont projetes sur 

10 la droite D2, etc. On a egalement represente sur la Fig. 5 une enveloppe £ de la 
constellation qui correspondrait a ce critere. Bien entendu, ce dernier se generalise h 
un reseau de dimension quelconque. Ainsi, pour une constellation parallelepipedique 
d'un reseau de dimension 3, le choix optimal correspondrait a une projection 
orthogonale sur une enveloppe polyedrique, le polyedre etant en fait un 

15 parallelepipede de faces paralleles a celles de constellation et dont tous les angles sont 
biseautes. 

Selon un second mode de realisation de Tinvention, on eflTectue la projection 
orthogonale du point re9U directement sur un sous-espace affine delimitant la 
constellation (par exemple sur Tune des droites Di,D2,D3,D4 ou Tun des points 

20 Pj2,P23,P34,P4i dans le cas de la Fig. 4) et Ton procede ensuite a une detection par 
spheres a I'interieur ce sous-espace afBne, Tavantage etant de travailler dans un 
espace de dimension N plus faible que la dimension k du reseau. 'Pour effectuer le 
decodage dans le reseau contenu dans le sous-espace affine, il faut projeter le point 
re9u directement sur le sous espace-aflBne. On pent neanmoins decider de n'eifectuer 

25 la projection que si le point re9u est sufiBsamment loin de la constellation. Les 
conditions de choix de projection sont done les memes que prec^demment. 

Nous montrerons tout d'abord que les points du reseau A de R*^ contenus dans 
le sous-espace afBne de projection, peuvent etre modelis6s par des points d'un reseau 
A' de . Une teUe mod^lisation permet d'effectuer le decodage par spheres dans un 

30 espace de dimension N<k et d'accelerer la detection. 

S*il existe L dimensions telles que pt^pT ou pt>^Pt , on fixe ces valeurs de pt 
a A/7 ou Mt selon les cas et Ton projette le point refu sur le sous-espace affine de 
dimension N-k-L correspondant. Sans perte de generalite, on peut faire Thypothese 




Sou z Vg la projection sur te sous-aspace aflfae de diaK=nd„„ du vecteur 
. D-ap.^ de p^JecUon .e> ,ue p„„ . ^ 'J^ 

-> done ecnre : ^ 

2''=bG+n=b''BVb''B+n 

(38) 

10 „nce de tai„e (^..)d.fi„ies co„,.e p..c.de„™e„t et „ est J^l: ' 
La simplification de la detection presentee dans cette section suppose aue le 

-l-5- ce7e^ 

vecteur liZt^ " " ^ '^^^-tion du 

vecteur hgne 2 appartenant au sous-espace vectoriel Tet nht^™. « 

vecteur constant z«=b°B : soustraction du 

2''=2''-2''=b''B+n 

20 (39) 



La matrice B contient vecteurs de base de taiUe k: z" Hn« 
bruite d'un reseau A' de dimension N de Pour ,ue la Ithod H 7 " """^ 

X , cspace K . 11 taut done trouver une matnVf. rp> 

25. contenant N vecteurs de base de taille N et p^n^^r^n. - 

genere par B dans Sait^t ^ '^'"^^ ^''"'^^'^"^ ^ 

^ P i» dans R . Sort B la matnce de taille (NxN) obtenue nar 

decomposition de Cholesky de la matrice de Gram BB^ Cette d 
garantitque: ' ^^"^ decomposition 

^0 BB'" =B'''B^''^ 

(40) 

■nfenenr^ Cette particular* pennet de feire r^conomje. de b. 
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decomposition de Cholesky normalement effectuee en debut de la d6tection par 
spheres. 

Soit U la matrice de passage de taille (tcxN) telle B'''=BU . 
Montrons que U = B'' (b'''"' J . En efFet , on obtient d'apres (40) : 

5 

BU=BB"(B'''-]r=B'''B'*'^(B"'-']r=B''' (41) 

Si on note x''=b''B le vecteur ligne de taille k representant un point du reseau 
A' dans R'^ , alors le vecteur ligne de taille N representant un point du reseau 
10 dans R" s'obtient par : 

x''=x''U=b''BU=b''B''' (42) 

On doit done feire subir 4 la meme transformation avant detection, h savoir : 

15 

z''U=(b''B+n)u=b''B'''+n' (43) 
La naatrice de covariance du vecteur de bruit n' est la suivante : 
20 4»''^n']=U'"4n''B]u=iVi)U'"U (44) 

soit encore, grSce k Texpression de U : 

4n'''n']=MB'''-'BB^(B'''-')r=iVoB'''-'B'''B^'^(B'''-')r=A^oIw (45) 

25 

Aucun blanchiment supplementaire du bruit n'est done necessaire apres projection. 

Les coordonnees entieres non fixees (les coordonnees fixees sont donnees par 
b") sont regroupees dans b'' et obtenues par une detection par spheres du vecteur 

30 Bien que I'invention ait ete expos^e dans le contexte d'un ^steme de 

telecommunication mobile DS-CDMA, elle s'applique egalement a un systeme de 
telecommunication mobOe MC-CDMA (Multi-Carrier Code Division Multiple 
Access). 
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Rappelons que la technique MC-CDMA combine la modulation OFDM 
{Orthogonal Frequency Division Multiplex) et la technique d'acces multiple CDMA 
Contrairement a la methode DS-CDMA dans laquelle le signal de chaque utilisateur 
est multiplie dans le domaine temporel pour etaler son spectre frequentiel, la signature 
multipUe ici le signal dans le domaine frequentiel, chaque element de la signature 
multipliant ie signal d'une sous-porteuse differente. 

A la reception les echantiUons du signal re9u sont soumis a une FFT et les 
signaux relatife aux differentes sous-porteuses sont ensuite filtres par un fihrage 
adapte a la signature de I'utilisateur et au canal de transmission. 

On a montre dans la demande fran^aise N° 0016351 deposee le 13 12 2000 par 
la demanderesse et incorporie ici par reference que si Ton note r(/) = (r,(/),...,r, (/)) 
le vecteur des signaux refus sur les L differentes sous-porteuses et d(/)=(^,(,%.,'2(0) 
le vecteur des K symboles 6tms par ou pour les ArdifKrents utilisateurs h I'LLt / 
on pouvait ecrire dans le cas oii les utilisateurs sont re9us de manidre synchrone : 

r(0=d(OAC(/)+T,(/) ^^g^ 



ou C(Oest une matrice traduisant I'etalement par les signatures des differents 

utdisateurs et I'effet des differents canaux de transmission, ti(0 = (^.(0 est 

20 un vecteur de bruit blanc additif gaussien et A est une matrice diagonale 
Diz^a^.^^) formfe par les anq>litudes des signaux des ou pour les differents 
utilisateurs. 

Comme indique dans la demande susmentiomiee, on peut montrer que la 
connaissance du vecteur d'observationy(/) = iy,(i),...,yj,(i)) ou : 

y(0=r(OC"(0 

Zc^^^sdcT ^'"'^ ^ ^" "^"^ vraisemblance du 

On en deduit que le vecteur d'observation peut s'expruner sous la foime: 

y(/)=d(/)AC(/)C''(i)+n(/)=d(i)M(0+n(/) (48^ 

oil B(i) est un vecteur de bruit colore et M(/)=C(OC"(0 
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Les vecteurs y(0, d(0, n(0 ainsi que la matrice M(0 figurant dans I'equation 
(48) sont a composantes complexes. La relation (48) peut encore s'6crire sous la 
forme reelle equivalente: 



y2(0=d2(/)M2(0+ii2(0 



(49) 



avec 



y20>iyr(i\y{(fX-^yl^(i\ylc(f)) on yJ^iflyiif) sont respectivement la partie r^eUe et la 
partie imaginaire du symbole yk(i); 

d2(0=(^,«(iV/(4 -.^K»M(0) di?(f\dl(i)sont respectivement la partie rtelle et la 
10 partie imaginaire du symbole <i*(/); 

m(/)=(«/'0'W(4 viJ?('W(0) oil n§Q\nl(f)sotA respectivement la partie r6elle et la 

partie imaginaire denaii) ; 

et ou Ml est la matrice IKxlK d^finie par : 



15 



-MJcM'—MkKM&c 



(50) 



avec Aftt=M5+ JMU on I'indice / a et6 omis pour simplifier les notations. 

L'equation (49) a une forme identique k celle de I'equation (6) et par 
consequent un algorithme de detection par spheres peut s'appliquer au vecteur 
d'observation, y compris la simplification proposee par I'invention lorsque le point 
re9u est situ6 en dehors de la constellation. 
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REVENDICATIONS 



1) Melhode de detection d'tme phuaiM Kie sy^boles (4(0) tt^mmis par o„ 
pour una ptalW AT d-utilisateu,. , partir d'un signal re^u. cha<p,e syn^ole d'un 
uctoteur appartenant 4 une comtelbtion de nodulation, fa ^t>od. de detection 

I^T, " * <^ 8'"*'' -naellations de 

.nod„la.K,„. ladtte plurality de symboles des differents utilisateu,. etant ,ep,^«e 

l»r ,^ po™ pamii u„ sous-ensemble de points dudit rfaeau, di, constellation et le 

^ re^u 4tant re^t* par un point caracteristi,ue de ce signal, di. point rc,^ 

^ d un p„„, de bdite constellation par un vecteur de bruit („). caracterfa«e ^ 

2Z J°'n "^"^ * l^"' - ™ -us- 

ypace affin^ d« de projecdon. paralBle ou confondu avec „„ sous-espace affine 

de^tant lad«e co„ste«on et une *ape de recherche du plus ptocj^ois^^ 
pomt amsi projet* paimi les points de ladite constellation. 

15 2) Methode de detection selon la revendication 1. caractWsfe en ce que led» 

sous-espace afflne d.Wta„, fa constellation es, d^tennin^ en fonction de la position 
du point re9u par rapport 4 ladite constellation. 

3) Methode de detection selon la revendication 2, ca,act«risfe en ce la 
^he du plus proche voisin es. Ii„*6e aux points de la constellation appatte,^ 
a une sphere centreesurledit point projett 

4) Mithode de dftection selon la revendication I ou 2. caract^risee en ce que Ic 

25 * ™ ---pace 

25 delunttan la constellation, la recherd^ du ph. proche voisin est effec«X paj^ 

pomtsdelaconstellationappartenantauditsous-espaceafflne. 

5) Mfthode de detection selon la revendication 4. caract^risfc «, ce que h 
recherche du plus proche voisin es, Wtee aux points du sous-espace affij 

0 aPI-enantatntesphSreduditsous-espacecentr^esurfeditpointprojete 
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6) Methode de detection selon Tune des revendications prec6dentes, 
caracterisee en ce que I'^tape de projection n'est effectuee que si le point re9u est 
eloign^ de ladite constellation de plus d'une distance pr^daerminee. 

5 7) Methode de detection selon I'une des revendications prec6dentes, 

caracterisee en ce que, les symboles de chaque utilisateur feisant robjet d'une 
multiplication par une signature de cet utilisateur avant d'etre transmis sur un canal 
de transmission, les coordonnees du point refu sont obtenues par une etape de 
filtrage adapts du signal re9u, le filtrage etant adapts aux canaux de transmission et 

1 0 aux signatures des dififerents utilisateurs. 

8) Dispositif de reception pour systeme de telecommunication DS-CDMA 
comprenant un dispositif de detection adapte a mettre en oeuvre la methode de 
detection selon I'une des revendications pr6cddentes. 
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9) Dispositif de reception pour systeme de telecommunication MC-CDMA 
comprenant un dispositif de detection adapte a mettre en oeuvre la methode de 
detection selon I'une des revendications 1 a 7. 
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